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Abstract 

We analyze the photon rates from a hadronic gas in equilibrium using chiral reduction formulas 
and a density expansion. The chiral reduction is carried to second order in the pion density which 
in principal includes all kinetic processes of the type X — > ttj and X — * nir^f. The resulting 
rates are encoded in the form of vacuum correlation functions which are amenable to experiment. 
The hadronic rates computed in this work along with the known perturbative QGP rates are 
integrated over the space-time evolution of a hydrodynamic model tuned to hadronic observables. 
The resulting yields are compared to the recent photon and low mass dilepton measurements at 
the SPS and RHIC. Predictions for the LHC are made. 
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I. INTRODUCTION 



Electromagnetic emissions in ultra-relativistic collisions are thermally dominated in the 
low to intermediate mass (M) and qx regions. Thermalization at RHIC has now been 
established from detailed flow measurements of hadrons. PHENIX has observed a large 
dielectron emission in the mass region below the p (M « 770 MeV). These emissions are 
much larger than those reached theoretically. The dielectron excess reported by PHENIX is 
one of the most dramatic piece of data stemming from RHIC. 

The PHENIX collaboration has also recently reported on photon spectra in the inter- 
mediate qx region. There is one caveat, in that PHENIX did not actually measure real 
photons but low mass dielectron spectra which was then extrapolated to the photon point. 
In this case their finding was an excess of photons at intermediate qr (1 < (?t GeV ^2.5) 
momentum. In this case the excess was in line with those reached theoretically. 
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Our qualitative summary of the above two PHENIX measurements is as follows. PHENIX 
has observed an excess of di-electron pairs at low mass. This excess is consistent with 
theoretical expectations at intermediate (1 < qr GeV < 2.5) momentum. The excess 
seen at low momentum < 1 GeV is well above any theoretical calculations done to date. 

Thermal emissions at low to intermediate invariant mass (M) and qx are involved due to 
the many reaction processes involving hadrons and the strong character of their interactions. 
The only organizational principles are broken chiral symmetry and gauge invariance, both 
of which are difficult to assert in reaction processes with hadrons in general. If hadrons 
thermalize with the pions and nucleons as the only strongly stable constituents, then there 
is a way to systematically organize the electromagnetic emissivities by expanding them not 
in terms of processes but rather in terms of final hadronic states. The emissivities are 
then amenable to spectral functions by chiral reduction. These spectral functions are either 
tractable from other experiments or amenable to resonance saturation. 

In section II, we derive the photon emission rates from a thermal hadronic environment 
in terms of one and two pion final states. We use the chiral reduction formulae to rewrite the 
rates in terms of spectral functions. In section IIA and IIB we compare our rates to some 
key processes based on kinetic theory as well as the leading order QGP rates. In section 
IIC we comment on the experimental extrapolation procedure used recently by RHIC to 
measure the photon emissivities by extrapolating the dilepton rates to the photon point. 

The remaining part of this works consists of comparison with experimental data. This is 
done by integrating the rates over the space-time evolution of a hydrodynamic simulation 
described in section III. Sections IIIA, B, C show our findings for WA98 and PHENIX. 
Section HID contains our predictions for the LHC. 

II. HADRONIC PHOTON RATES 

For a hadronic gas in thermal equilibrium the number of photons produced per unit four 
volume and unit three momentum can be related to the electromagnetic current-current 
correlation function hi 




a, 



471-2 



em 
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with q 2 = and 

W(q) = J d 4 x e- iq - x Tr ( e - {H - F)/T 3 fl (x) J M (0)) . (2) 

In the above expression J M is the hadronic part of the electromagnetic current, H is the 
hadronic Hamiltonian and F is the free energy. Below the phase transition the trace is 
carried over stable states with respect to the strong interaction {e.g. pions and nucleons). 
From the spectral representation and symmetry we can re-express the correlator in terms of 
the absorptive part of the time-ordered correlation function 

W(9) = TT ^ImW F ( g ), (3) 



where 



W F (q) = i j d 4 x e iq - x Tr (e~ {H - F)/T T3^(x) J„(0)) . (4) 

In this work we will consider a heat bath which is nucleon free. In this case the trace can 
be expanded as 

W F (q) = W + j d^W, + 1 J dTTidTTaW™ + • • • , (5) 

where di\i are the pion phase space factors given by 

d 3 ki n(Ei) 

In the above density expansion we have defined 

W = i f d 4 xe tq - x (0\TJ^x)3^0)\0) 



W, = % j d 4 xe iq - x (TT a {k 1 )\T3^x)3 fl {0)\7r a {k l )) 

W™ = ij d 4 xe* q - x (K a {k l )n\k 2 )\T3^x)3^ a (h)*\k 2 )) (7) 

where the latin indices a, b are summed over isospin. These are the first three terms in an 
expansion in terms of the pion density 1 . 



1 More specifically the dimensionless expansion parameter is k w n v / (2m ff /^) where n, is the pion density. 
This corresponds to k w 0.18,0.30,0.84 at temperatures of 120,140 and 190 MeV. We therefore expect 
the expansion to be reasonable unless new thresholds open up. 
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The first contribution in ((Tj) is dominated by the transverse part of the isovector correlator 
and is fixed entirely by the measured electroproduction data. It vanishes for real photons 
since the heat bath is stable against spontaneous photon emission 2 . Therefore W = 
and does not contribute to the photon emissivities. The next two terms, W,,. and ~W nn , 
can be reduced to measurable vacuum correlators by the chiral reduction formulae [2j. The 
one-pion reduced contribution is by now standard. Its dominant contribution involves 
VV and AA correlators in the vacuum and reads |3| 

12 

F 

J 7T 



W%(q,k) = -q 2 lmn v (q 2 ) 



^(k + q) 2 lmU A ((* + q) 2 ) + (q -> -q) 

J TT 



+ — {{k ■ qf - mlq 2 ) ImIIy(g 2 ) x ReA R {k + q) + (q -> -q) 

J TT 



(8) 



where ReA# is the real part of the retarded pion propagator, and Ily and II^ are the 

transverse parts of the VV and AA correlators. Their spectral functions are related to both 

fl 

e + e annihilation and r-decay data as was compiled in [4|. 

The two-pion reduced contribution is more involved. Its full unwinding can be found 
in We quote only the dominant contributions 
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hv^{q,k x ,k 2 ) = (2k 1 + q)„k lv ReA R (k 1 + q)}Im%% (q,k 

J TV J TV 

+ (q ~q) + (h -> -ki) + (q, h -> -q, -ki) 

+ -1^(2^1 + qYReAnih + q)e a3e e^lmB;i(k 1 , k 2 ) 

J TT 

- 1 [<T - (h + qr{2h + qTReAnih + q)} 

J TT 

x e «3e e a3/ Imj ge/^ i + q ^ 2 ) 

+ + qTih + q y(2k 1 + q) 2 [ReA R {k x + q)} 2 

J TT 

x e a3e e a3f lmB e fl i(k 1 + q,k 2 ) + (k 1 ^ -h) (9) 



TT'y J 



The first term contains the pion-spin averaged 7T7 forward scattering amplitude (Im7^ 
For real photons (q 2 = 0) this is entirely constrained by photon fusion data by crossing 
6| while for virtual photons it can be worked out by further chiral reduction as 



symmetry 



It also vanishes for massive photons having M < Im-n- 
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FIG. 1: Pictorial chiral reduction of B. 

done in The remaining terms all contain the contribution labeled B which reads 

B;t(ki, k 2 ) = ij d 4 xe^ x (7rUh)\T (j^(^)ji(O)) \7ti(k 2 )) , (10) 

and is still amenable to further chiral reduction. The dominant contribution arises from 
the schematic diagrams shown in Fig. [1] resulting in the following spectral contributions 

ImBjKh, k 2 ) = A**/ [g^ kl + k 2 f - {h + kd^h + fc 2 )„] Imn y ((h + k 2 f) 

J TT 

+ (h -> -k 2 ) - ^6 ef [g^kj - k lfl k lu ] lmU A (kj) (11) 

J TT 

It is also instructive to look at the above results in the limit where the incoming pions are 
soft (ki, k 2 — > 0). We then find the simple result that 

1 1 fi 

74 W -(«) = 71 ( Imn A(g 2 ) - ImlW)) (12) 

J TT J TT 

which reproduces the result of Eletsky and Ioffe|7]. We should also mention on a more 
practical level that the dominant contribution to the above expression comes from terms 
containing II y and one can therefore neglect terms in Eq. ([9]) containing 7^ 7 or ILi. The 
dominant mechanism is therefore the first diagram on the right hand side of Fig. [Q This 
corresponds to the matrix element for tttt scattering via exchange of a vector meson. The 
additional terms in Eq. (Q attach a photon to one of the external lines. 

In Fig. [2] we show the photon rates stemming from and for three temperatures. 
Clearly the 2-pion emission rates provide substantial enhancement when the emitted photon 
is soft due to pion bremsstrahlung. At the highest temperatures, ~W n7r provides an enhance- 
ment even for high energy photons. But at these temperatures the virial expansion is clearly 
beyond its limit of applicability. We stress that the hadronic contributions to involves 
all hadronic processes in the heat bath with two pions and a photon in the final state. 
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FIG. 2: Photon emission rates and W nw for T = 100, 150, 200 MeV. 
A. Comparison to Kinetic Theory 

One of the nice features of the above spectral function approach is that one does not need 
to organize the calculation in terms of the many possible reaction mechanisms, as there is 
no a priori small expansion parameter in the strongly coupled hadronic phase. Instead the 
calculation is organized in a virial-like expansion and all possible reaction channels should 
be included via the zero temperature spectral densities. The drawback is that we loose some 
physical intuition as we tend to discuss photon production from a kinetic theory approach. 
In this section we now compare our results with kinetic theory results performed by a number 
of other groups. 

n 

First let us start with our W„- piece, which was originally worked out in |3J for both 
dileptons and photons. From a kinetic theory standpoint this should include all reactions 
that contain one pion in the final state. An example of the dominant contribution in this 
channel is 7rp — > -k^. In Fig. [3] we show the one pion piece W,,- as well as the rates stemming 
from the reaction up — > tt^ computed in |8j. On a qualitative level we find very good 
agreement considering the rates were derived from two very different approaches. We should 
mention that a full SU(3) calculation [9] also includes reactions containing K and K* (e.g. 
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FIG. 3: Photon rates from spectral functions versus kinetic rates. See text for details. 



nK* -> AT 7 ). 

We now come to the piece which is new to this work. From a kinetic theory 

standpoint this should be identified as including all reactions with two pions in the final state 
and in principal can contain any number of initial state particles. We expect the reactions 
7T7T — ► p7 and p — > txtx^ to dominate. These reactions were calculated in [8|] and are shown in 
Fig. [3j These two kinetic processes alone cannot explain the observed enhancement that we 
find in our contribution. Our computed ~W nn should also contain 2 — > 3 processes such 
as 7T7T — > 77777. This reaction has been worked out beyond the soft photon approximation 
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10|. We show a parameterization of their results in the region 0.1 < go (GeV) < 0.5 in 
Fig. [3j On a qualitative level we see that the pion bremsstrahlung agrees with our 
contribution. 

We do not claim that the rates encoded in W„- and W W7r are the end of the story. In the 
above comparison we have selected the reaction channels which we expect to be implicitly 
included in our analysis. We can conclude that we have a good handle on many of the 
kinetic theory processes shown above. The most significant being the 7777 bremsstrahlung 
encoded in W^. There are other reactions (e.g. uj — > 777), while not included in our SU(2) 
analysis, have been shown to have strength at low energy 



S 



B. Comparison to QGP Emission 



There has been great progress in the calculation of the QGP photon rates in perturbative 
QCD. We will not go through the full history of these results but instead highlight some key 
points as they may become relevant from a phenomenological perspective. 

First off, the born contribution corresponding to qq — > 7 contributes at order a® for dilep- 
tons but vanishes at the photon point due to energy momentum conservation. Therefore, for 
photons, the first non-trivial contribution will come from one-loop diagrams, 0(g 2 ), corre- 
sponding to the annihilation q + q — > (77 and compton g + q(q) — > q(q) +7 processes. These 
diagrams contain a logarithmic singularity stemming form the exchange of a soft massless 
quark. The singularity is cured by including HTL 12J corrections to the quark propagator. 
The net result is to give the exchanged quark a thermal mass of order gT and the rates are 
indeed finite as shown in ll|, [l3]. We have plotted the photon rate due to 2 — > 2 process in 

"V- u. — ,. HH B — . — — — <»— 

scattering are not the only diagrams that contribute at order g 2 . It was found that there 
are a class of diagrams that even though they are naively of higher order are promoted to 
order g 2 due to co-linear singularities. Two examples of such processes are bremsstrahlung 
(qx — > jqx) and annihilation with scattering (qqx — > jgx) where x can be either a quark, 
ant i- quark or gluon. 

In order to obtain the leading order result one must resum an infinite set of ladder 
diagrams which contribute at the same order. In resumming this infinite set of diagrams 
one must also take into account the LPM effect which suppresses the rate due to the effects of 
multiple collisions on the photon production process. This was demonstrated and evaluated 
in full in , 119] . In this work we use the complete leading order rates including the 



LPM effect as parameterized in [18J. The additional processes included in the full leading 
order calculation give a large enhancement to the rates as seen in Fig. [5j 

The leading order QGP rates have one parameter, a s , which should be evaluated at a 
scale on the order of the temperature. In this work we evaluate the coupling constant using 



the two-loop (3 function in MS scheme. In order to asses the uncertainty in the scale choice 
we have evaluated the yields assuming the coupling runs with two different scales: \x = tiT 
and \i = 2ttT. At T = 400 MeV, which is a typical initial temperature in our hydrodynamic 
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evolution we have a s (7iT) = 0.38 versus a s (27rT) = 0.26. While at our chosen transition 
temperature of T = 190 MeV we have a s (7rT) = 0.75 versus a s (27iT) = 0.40. 




FIG. 4: Scale dependence of the perturbative QGP yields at SPS, RHIC and the LHC. 

Even though we have yet to discuss the details of the space-time evolution (see Section HTT|) 
it is still instructive at this point to compare the yields from the QGP phase for various 
systems in order to see the uncertainty in the scale choice. In Fig. H] we show the ratio of the 
QGP yield coming from collisions at SPS, RHIC and the LHC at the two scale choices. We 
will focus on the results above q± ~ 1 GeV since the perturbative rate calculations break- 
down below this scale as we will discuss in a moment. In general we find an enhancement 
in the yields by a factor of ~ 60% at the SPS and ~ 30% at the LHC by decreasing the 
scale from 2irT to tcT. At the LHC there is a smaller uncertainty in the scale choice since 
the temperatures in the system are much higher yielding smaller running coupling effects. 
This analysis shows that at RHIC there will be a systematic uncertainty in the yields on the 
order of ~ 40% due to the scale choice. With this in mind we have chosen to fix the scale 
to fi = 7rT for the remainder of this work in order to get an upper bound on the photon 
emissivities from the QGP phase. 

Let us now come to Fig. [5] where we have compared the Hadronic and quark-gluon plasma 
rates for a broad range of photon energies. We have generated the plot at the transition 
temperature T = 190 MeV, where both rates have the largest theoretical uncertainties. The 
QGP rates are plotted with a substantial coupling a s = 0.75 which is the highest reached 
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FIG. 5: Hadronic rates from the spectral function approach versus the perturbative QGP rates. 



in the simulation. We should stress that in the transition region one should be extremely 
skeptical of the rates. For both the hadronic gas and QGP the expansion parameters have 
become of order one (k ~ 0.8 and g ~ 3). One might argue that the QGP rates can be 
trusted as long as the argument under the log is larger than one. This is the case as long as 
one remains above go ~ 1 GeV. It is reassuring to see that the photon rates are in qualitative 
agreement (within a factor rs 2) as one approaches the transition region from above or below. 



C. Photons versus low mass dileptons 



In this section we would like to discuss some general considerations between photons and 
low mass dileptons as well as the rates for the latter. The discussion of this section will be 
relevant for understanding the systematic uncertainty in the recent photon measurements 
at PHENIX as well as understanding part of the source of the low mass di-electron excess 
observed at PHENIX. We will withhold the model comparison to data until Section IIII CI 

The rate of dilepton production per unit four volume and four momenta is given by 

dR -a 2 ( 2m, 2 \ ( 4m 2 \ 1/2 1 p 



n 



for massive dielectrons. As a result we have 



dR 2a ( 2m? \ ( Am 2 \ 1/2 / n dN 



#, ^\ 1+ f)v—f) l 9 wJ (14) 

which ties the dielectron rate to the virtual photon rate N* for spacelike momenta. In the 
above expression the correlation function W(g) is the same as in Section [TT]but now evaluated 
at q 2 < 0. Therefore the spectral function approach used for the rates in the hadronic phase 
naturally contains the rates for photons and dileptons. The resulting dilepton rates when 
including the additional two pion piece can be found in jsj and will be briefly discussed 
later. 

The recent photon measurement reported by PHENIX is actually extracted from low 
mass dielectron data below the two pion threshold. This is done in order to avoid the 
large background from hadronic decays. In order to extract the photon rates the PHENIX 
collaboration makes use of the following relation between photon and dilepton production 

In the above expression S is a process dependent factor. This process dependent factor goes 
to 1 for q 2 —>■ 0. For 7r and rj decays (which is the dominant hadronic background) S — > 
when M > M/^the mass of the hadron). By cutting out invariant masses less than the it 
and rj mass the background from hadronic decays will be largely suppressed. 

The assumption taken by the PHENIX experiment is that S ~ 1 when M <C q?- The 
high q?, low mass electron pairs are then taken to the photon point via Eq. (fT5l) . Since 
our spectral function approach has direct access to both photons and low mass dileptons 
we can easily determine (in a model independent way) the S factor for hadronic production 
processes. 

In Fig. [6] we plot (in arbitrary units) the virtual photon rate (q°dN* /d 3 q = Sq°dN/d 3 q) 
where q°dN* /d 3 q is the rate of virtual photon production via Eq. f|T4l) . The left figure 
includes only W,,- while the right figure includes W,,- + W^. We have fixed the three 
momentum to q = 0.25, 0.5, 1, 3 GeV. If our process dependent factor S was indeed equal 
to one, we would not see any mass dependence in our plotted q°dN* /d 3 q. This is not the 
case however. For our leading order (W,,-) rates at q > 0.5 the factor S varies by a few 
percent up to a mass of 300 MeV. If this was the dominant thermal process then the photon 
extrapolation made by PHENIX would be satisfactory. However, when including the more 
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complicated production mechanisms contained in W T7r the S factor can change by as much 
as 30% at M — 300 MeV when q — 1 GeV. Since PHENIX fits the low mass electrons in the 
100 < M < 300 MeV region we can confidently say that our S factor will give a discrepancy 
of « 15%. This leads to a fairly large unconstrained systematic error on the PHENIX 
extrapolation procedure. We stress the fact that this error is unconstrained because the S 
factor could in principal be very different in the QGP where the production mechanisms 
may be much different. 
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FIG. 6: Rate of virtual photons (arbitrary units) from the spectral function approach including 
(left) and + W W7T (right). 



Let us conclude this section by a quick discussion on the dilepton rates in the QGP 
phase. These will be used in our analysis of the PHENIX data in Section fill Bl Naively one 
would expect the leading order Born contribution, qq — > 7*, to dominate the dilepton mass 
spectrum. This breaks down at low mass due to phase space constraints. For M < m t h 
the next to leading order contribution starts to take over. All of the same considerations 
for photons has to be taken into account for dileptons including the HTL resummation, 
the summation of ladder diagrams and the inclusion of the LPM effect. This was worked 
out in full for dileptons in [20j. For applications we use the parameterization of this result 



by 2l|, |22j. These rates are applicable in the kinematic regions where go > M, g > T 
and M 2 ~ g 2 T 2 . The PHENIX acceptance can fall outside of these regions, so we caution 
the reader to this additional systematic uncertainty that comes about when comparin g to 
experiment. To the LPM rates one must also add the 2 —>■ 2 process worked out in 23] 
which were computed in the approximation M < T and go ^ T. 
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III. SPECTRA IN ULTRA- RELATIVISTIC HEAVY-ION COLLISIONS 



Before a direct comparison with data can be made the photon emission rates of the 
QGP and Hadronic phases must be integrated over the space-time evolution of the collision. 
The collision region is modeled using a relativistic hydrodynamic simulation tuned in order 
to reproduce hadronic observables. In this section we discuss the specifics of the model, 
including the initial conditions and equation of state (EoS), but leave the technical details 
to the literature. 

The required initial condition for the hydrodynamic evolution is the entropy density, 
s(x,y), and baryon density, ng(x,y), in the transverse plane for a given impact parameter 3 
(b) and initial proper time (r ). The entropy and baryon density are proportional to the 
number of wounded nucleons (participants) which is set by a Glauber model. More precisely 
the initial condition is 

7~o CL T 

i u\ C B dN WN 

n B (r,b) = — (r;6), (16) 

To a r 

where C s = l/N part dS/dy (Cb = ^/N part dn B /dy) is the entropy (net Baryon) density per 
unit rapidity per participant and is set in order to reproduce the charged particle multiplicity 
and net proton number respectively. Once the initial conditions are set the hydrodynamic 
evolution equations are solved. The resulting evolution yields the energy density and flow 
velocity as a function of proper time r and transverse coordinate. The energy density can 
then be converted to a temperature using the given EoS. The rates are then integrated over 
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25] 



the space-time volume following the same procedure in 

Table [I] shows the parameters used for the SPS, RHIC and LHC. For RHIC we have used 
two different evolution models (labeled as RHIC 1 and 2) in order to study the effects of the 
space-time evolution. We will discuss the specific evolution models as they become relevant 
in the following sections. 

We should stress that we make no attempt to perform global fits to the data. There are 
still a number of uncertainties present in the hydrodynamic model which would effect the 



In the present work we actually do not perform off-central calculations. Instead the simulations are done 
in 1+1D where the system size is fixed in order to have the same number or participants. The resulting 
yields agree on the 10% level and are therefore within the overall uncertainty in our hydrodynamic model. 
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Parameter 



SPS RHIC 1 RHIC 2 LHC 



y/a NN [A-GeV] 


17.3 


200 


200 


5500 


A 


208 


197 


197 


208 


^NN [ mb ] 


33 


40 


40 


60 


C s 


8.06 


20.8 


20.8 


42 


C B 


0.191 


0. 


0. 


0. 


EoS 


BM 


Lat 


Lat 


Lat 


Centrality: 


0-10% 


0-20% 


0-20% 


0-20% 


b [fin] 


3 


4.5 


4.5 


4.8 


Npart 


340 


269 


269 


293 


t [fm/c] 


1 


1 


0.5 


0.5 


T(r_L = 0, r ) [MeV] 


245 


336 


398 


501 


Tfrzout [MeV] 


120 


140 


160 


140 



TABLE I: Hydro dynamical parameters for: SPS, RHIC and LHC. 

resulting photon yields. The approach taken in this work is to choose parameters for our 
evolution model which yields a reasonable description of bulk observables such as the total 
multiplicity, pt spectra and elliptic flow. An example of one uncertainty is the shear viscosity, 
which when included would require us to re-tune the initial conditions in order to achieve 



the correct final state multiplicities (as shown in 2a]). In addition, viscous corrections will 



also modify the underlying photon and dilepton rates [27j, [28j. These considerations are all 
beyond the scope of this work and should be thought of as part of the overall uncertainty 
in our evolution model. 



A. SPS: photons at WA98 



3+Pb collisions at yfs = 158 
I30I ]. In this case upper 



The WA98 collaboration has measured direct photons from P 
GeV per nucleon using the statistical subtraction method in 
limits could be obtained for photon momentum in the range 0.5 < q±_ < 1.5 and data 
points obtained above 1.5 GeV. In addition the yield of direct photons was also measured 
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3l| using direct photon interferometry. The most probable yield is found by assuming a 
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source size of 6 fm and a lower limit is obtained by assuming a vanishing R ou t- The data is 
summarized in the left plot of Fig. [7J 

Fig. [7J shows our predicted yields using the SPS hydrodynamic evolution from Table fl] 



In this case we have used identical parameters to [32]. This work used a bag model EoS 
having a latent heat of 0.8 GeV/fm 3 and was able to reproduce the p? spectra and elliptic 
flow at the SPS. Above 2 GeV one must also include the prompt production which we have 

n 

left out in this analysis. This additional contribution was discussed thoroughly in [8J. Our 
results are consistent with the upper limits that were obtained at intermediate momentum. 




FIG. 7: Photon Spectra at the SPS. (Left) Data comparison for < q± (GeV) < 4. (Right) Same 
as left but for < q± (GeV) < 0.8. See text for details. 

The right figure of[7Jis the same as the left but rescaled in q T between and 0.8 GeV. 
This figure shows the leading order QGP along with the single pion process are almost 
an order of magnitude below the data. By including the additional ~W n7r component our 
rates are now at the lower limit of the data. In addition one may question how changes 
to the evolution might effect the results. It is very probable that one could find a different 
set of parameters which reproduces the hadronic data at the SPS while also increasing the 
yields of direct photons. This type of global analysis is beyond the scope of this work. It 
is clear however, that pion bremsstrahlung as encoded in is a necessary component to 
understanding the direct photons at SPS. 
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We should stress that the role played by baryons is large at low momentum [33|]. While 



baryons were not included in this work, we can sti 
other calculations. The work of Steele et al. 



34 



1 estimate the enhancement by looking at 



351 ] found that baryons, when treated to first 



order in nucleon density, can increase the rates relative to W„- by an order of magnitude at 
~ 0.5 GeV. Based on Fig. [2] we can estimate that the inclusion of nucleons could increase 
our overall rates (when the irN contribution is added to W^) by a factor of about two. 
This would comfortably keep us within the WA98 data. 



B. RHIC: photons at PHENIX 

Let us now come to the recent PHENIX measurements of photons. For Au+Au collisions 
at RHIC energies = 200 GeV) we have chosen to use two different evolution models 
(called RHIC 1 and 2) in order understand some of our systematic uncertainty stemming 



from the hydrodyanmic evolution. In both cases we use a lattice motived EoS [36] which 
includes a rapid crossover followed by an interpolation into the hadronic resonance gas phase. 
Even though there is no true phase transition, we still must choose a value of T c where we 
switch from QGP to hadronic production. Of course, this choice of T c will not effect the 
hydrodynamic evolution which is determined only by the sound speed. We will fix this to 
T r — 190 MeV which is consistent with lattice calculations at almost physical quark masses 



[37|, 



The difference between RHIC 1 and 2 is the initial hydrodynamic starting time as well as 
the freezeout temperature. For RHIC 1 we have r = 1 fm/c and Tf rzout = 140 MeV. In the 
case of RHIC 2 we have started the hydrodynamic evolution even earlier at tq = 0.5 fm/c. In 
this case, in order to have reasonable agreement with bulk observables one must use a higher 
freezeout temperature, Tf rzout = 160 MeV. The RHIC 1 parameter set is closer to the one 



used in 



26 



39( | which found good agreement with the elliptic flow of hadrons. In the case of 
RHIC 2 we expect to have a larger contribution coming from the QGP phase and a smaller 
contribution from the Hadronic phase. Note that we do not include any out-of-equilibrium 
contribution (other than the prompt contribution) before our initial starting time r . 

In Fig. [8] we show the hydrodynamically evolved photon emissivities using the evolution 
parameters of Table [I] for RHIC 1 (left) and RHIC 2 (right). The yields from the hadronic 
gas and QGP cross at q± w 1.5 and 0.6 for RHIC 1 and RHIC 2 respectively. The earlier 
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FIG. 8: Photon spectra from hadronic sources versus the QGP at RHIC. Left is for the evolution 
model RHIC 1 and the right is for RHIC 2 (see Table H]). 

thermalization in the RHIC 2 evolution gives a much larger QGP contribution at high 
momentum. 

We now come to the recent measurement of direct photons by PHENIX 40| . We will now 



be slightly hypocritical. In Section Ul Cl we have stressed that there will be a large uncertainty 
in extrapolating the photon signal from low mass dilepton data. We have demonstrated 
that for the hadronic gas used in this work that the discrepancy can be as large as ~ 15%. 
Regardless we will plot our results for photon production on-top of the PHENIX data keeping 
in mind that there is at least a 15% uncertainty in their extraction. 

Fig. [9] shows the evolved photon rates versus the RHIC data for the two different hydro- 
dynamical set ups in Table [J (left) is RHIC 1 and (right) is RHIC 2. We have now also 
included the prompt photon production which is derived from a fit to the measured photons 
in pp collisions scaled by the number of binary collisions. 

Our finding is that even though the relative contributions from the different phases are 
much different the sum remains about the same for the two different evolution models. 
In the case of RHIC 1 the relative contributions from the QGP and hadronic are about 
equal. For RHIC 2 the QGP has a larger contribution to the overall yield (which grows 
with increasing q±). The differences in the two models are largest at high q± where the 
prompt production dominates. Taking into consideration all of the systematic uncertainties 
that we have discussed throughout this work it is impossible to discriminate between the 
two evolution models. This suggests that there may be a large ambiguity in using photons 
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FIG. 9: Photon spectra at RHIC compared to the recent PHENIX data. The left plot is the 
evolution set RHIC 1 and the right is for RHIC 2 (see Tabled]). 

as an early time probe of the medium. 



C. RHIC: low mass dileptons at PHENIX 



We would now like to revisit our past analysis 41] of the measurement of dielectron 



pairs by PHENIX (42 



431 ] . In our prior work we used the leading order Born contribution 
(qq ^7*) as the only reaction present in the QGP phase. While this is true at high mass 
the naive perturbative expansion breaks down at low mass which might explain the missing 



low mass yield from theoretical models [44|. We have therefore included in addition to the 
Born contribution the next to leading order contribution. The rates used in this analysis 
were summarized at the end of Section III CI For the analysis of the low mass dileptons we 
have chosen to use the RHIC 2 evolution model. Even though this parameter set is for more 
central collisions b = 4.5 fm we will make a direct comparison with the min. bias data. In 
principal we could perform runs at various centralities and average accordingly in order to 
make a more direct comparison but this is beyond the scope of this work. 

In the left plot of Fig. [TU]we show the individual contributions to the low mass dielectron 
yield. New to this figure is the inclusion of the W W7r component in the hadronic phase. Even 
though this component added a lot of strength to the low momentum photons the effect is 
much less dramatic in the case of dileptons. There is a large enhancement in the dilepton 



rates at low mass and low qx as demonstrated in [5|. By M ~ 100 MeV almost all of the 
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FIG. 10: Dielectron spectra as measured by PHENIX. The left shows various individual contribu- 
tions to the thermal yield. The right figure shows the total yield using 1) the leading order born 
contribution and 2) the NLO QGP contribution. 



enhancement contained in is removed by the low My cuts of the PHENIX acceptance. 
Therefore, to a very good approximation, our W„- makes up the entire hadronic yield at 
PHENIX after the acceptance cuts. Also shown in Fig. [TU] is the qgp Born contribution. 
As we can see this contribution diminishes at low mass due to phase space. However, when 
we include the next to leading order (NLO) corrections, the QGP yields starts to increase 
below 2Mth an d actually overtakes our hadronic contribution below M ~ 500 MeV. At high 
mass (above ph 1 GeV) the NLO rates are suppressed relative to the Born term due to the 
thermal quark mass 4 . 

The right plot of Fig. [TUJ shows the total yields which includes summing the cocktail 
provided by PHENIX along with our thermal hadronic and qgp yields. We have shown the 
net yield using only the Born contribution ("Sum w/ Born only") which is essentially the 
result from our previous work [4l|. We also show the yield using the full leading order in 
a s result ("Sum w/ NLO QGP"). The leading order QGP explains the excess in the mass 
region 100 < M (MeV) < 150. About 50% of the thermal qgp production in this mass 
region is due to the one-loop Compton and annihilation processes. The other 50% required 
to explain the excess comes from the higher-loop processes included in the ladder diagram 



4 One could reproduce the NLO rates above ss 1 GeV by simply using the Born rate with massive (m ~ gT) 
quarks. 
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resummation. The two-loop piece of this corresponds to the bremsstrahlung diagrams. One 
should not be surprised that the QGP contribution overtakes the hadronic yield at low mass. 
It is already known that quark bremsstrahlung dominates over ir bremsstrahlung in the q 



integrated rates 



46j. In addition the qx spectrum of the tt bremsstrahlung falls much 
quicker then that of the quark bremsstrahlung , which only increases the relative yield of 
the quark to pion contribution when low qr cuts are applied. Larger transverse flow in the 
hadronic gas alleviates this picture by shuffling strength from low qj- to high qx but it can 
not compensate for the differences in the qr spectrum. 

Even though the addition of the full leading order qgp rates does not explain the excess 
in the mass region 250 < M (MeV) < 550 it does help in adding some strength to the 
yield in this mass region. The mechanism behind the excess in the 250 — 550 MeV region 
is still unclear. What the authors find striking is the resemblance between the excess seen 
here and the excess observed in lower energy C + C and Ca + Ca collisions referred to 
as the "DLS puzzle". The resolution was related to a better theoretical understanding of 
bremsstrahlung production in p + n collisions and verified by looking at p + p and p + d 
control measurements [4^]. In light of this, it might be useful to have a cocktail available 
which has been constrained by measured dileptons in p + d or d + Au collisions at RHIC 
energies. 



D. LHC: predictions 



In this final section we would like to make some predictions for the photon production 
rate we expect for LHC conditions. It is very difficult to tune our hydrodynamic model 
without having any data on the net multiplicities. Instead we have set our inital condition 
in order to approximately yield the predicted multiplicity 48|, |49(. In going from RHIC 
to LHC the inelastic nucleon-nucleon cross section, which enters into the wounded nucleon 
profiles has increased from 40 mb to 60 mb. We have also doubled the initial entropy density. 
Our final parameter set is summarized in Table [J and is in qualitative agreement with the 



hydrodynamic models used in 



50 



5l|. 



In Fig. [TT] we display our projected photon emissivities at LHC. The hadronic and qgp 
yields cross around qx ~ 1 GeV with the hadronic gas dominating at low momentum. We 
have also included predictions for the prompt photon production in p+p at LHC energies 
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FIG. 11: Prediction for photon yields from Pb+Pb collisions at y/s = 5.5 TeV. Left: Relative con- 
tribution from W,,-, Wtt,,- and the QGP phase for < qr (GeV) < 2. Right: Relative contributions 
from Hadronic (W,,- + W^), QGP and prompt production for I < qr (GeV) < 5. 



52j scaled to Pb+Pb collisions. The prompt production starts to dominate the yield above 



qr ~ 2.5 — 3.0 GeV similar to what was found at RHIC. 



IV. CONCLUSIONS 



In this work we have provided an analysis of the photon emissivities over a large range of 
collision energies ranging from the SPS (-y/i = 158 GeV) to RHIC (y/s = 200 GeV) and to 
predictions at the LHC (-y/s = 1.5 TeV). The analysis incorporates both hadronic and QGP 
rates integrated over an underlying hydrodynamic evolution tuned to bulk observables. In 
this section we would like to summarize the main conclusions of our work. 

• We have extended the chiral reduction approach to include terms to second order 
in pion density. The new addition of this work is inclusion of processes of the type 
X — > TTiT'y. The use of a spectral function approach allows us to treat both dilepton 
and photon production on an equal footing and incorporate broken chiral symmetry 
in a systematic fashion. We stress that our hadronic rates are parameter free and are 
completely constrained by r decay, electroproduction and Compton scattering data. 

• Using our hadronic rates, we have tested the low mass extrapolation used by PHENIX 
and found that we may expect a systematic error in their procedure of at least 15%. 
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• We have found that the pion bremsstrahlung process, encoded in W^, leads to a 
considerable photon enhancement in the low energy region. These processes are a 
necessary ingredient in a quantitative understanding of the low energy (q± < 500 MeV) 
photon data as demonstrated by the recent WA98 measurements (see Fig. [7]). 

• The hadronic and LO qgp photon rates are able to adequately describe the RHIC 
data. We have found that two different evolution models, which have different relative 
yields from the QGP and hadronic phases, are both able to describe the data. 

• We have re-addressed the low mass dilepton data from PHENIX. Our findings are 1) 
the W nn processes are completely removed by the acceptance cuts. 2) Inclusion of the 
NLO dilepton rates from the QGP are able to account for the enhancement found in 
the M = 100 - 150 MeV region. 
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